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Abstract. For any pair of quantum states, an initial state |/) and a final quantum state \F), 
in a Hilbert space, there are many Hamiltonians H under which |/) evolves into \F). Let us 
impose the constraint that the difference between the largest and smallest eigenvalues of H, 
£'max and i?min, is hcld fixed. We can then determine the Hamiltonian H that satisfies this 
constraint and achieves the transformation from the initial state to the final state in the 
least possible time r. For Hermitian Hamiltonians, r has a nonzero lower bound. However, 
among non-Hermitian T'T-symmetric Hamiltonians satisfying the same energy constraint, 
T can be made arbitrarily small without violating the time-energy uncertainty principle. The 
minimum value of t can be made arbitrarily small because for T^T-symmetric Hamiltonians 
the path from the vector |/) to the vector \F), as measured using the Hilbert-space metric 
appropriate for this theory, can be made arbitrarily short. The mechanism described here 
is similar to that in general relativity in which the distance between two space-time points 
can be made small if they are connected by a wormhole. This result may have applications 
in quantum computing. 

Key words: brachistochrone; FT quantum mechanics; parity; time reversal; time evolution; 
unitarity 
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1 Classical brachistochrone problem 

Three hundred years ago the solution to a famous problem in classical mechanics known as the 
brachistochrone was found almost simultaneously by a number of distinguished mathematicians 
including Newton, Bernoulli, Leibniz, and L'Hopital. (The term brachistochrone is derived 
from Greek and means shortest time.) The problem is stated as follows: A bead slides down 
a frictionless wire from one given point to another in a homogeneous gravitational field. What is 
the shape of the wire connecting the two points that minimizes the time of descent of the bead? 
The solution found by these mathematicians is that the wire must be in the shape of a cycloid. 

Of course, it is implicitly assumed in the derivation of the brachistochrone that the path of 
shortest time of descent is real. It is interesting that if one allows for the possibility of complex 
paths of motion, one can achieve an even shorter time of flight. 

To illustrate how shorter times can be achieved by means of complex paths, let us consider 
the simple classical harmonic oscillator, whose Hamiltonian is given by 

H = p^ + x^. 

If we have a particle of energy E = 1, then the classical turning points of the motion of the 
particle are located at x = ±1. The particle undergoes simple harmonic motion in which it 

*This paper is a contribution to the Proceedings of the Seventh International Conference "Symmetry in 
Nonlinear Mathematical Physics" (June 24-30, 2007, Kyiv, Ukraine). The full collection is available at 
|http://www.emis.de/journals/SIGMA/symmetry2007.html| 
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oscillates sinusoidally between these two turning points. This periodic motion is indicated in 
Fig. [T] by a solid line connecting the turning points. However, in addition to this oscillatory 
motion on the real-x axis, there are an infinite number of other trajectories that a particle of 
energy E can have [1]. These classical trajectories, which are also shown in Fig.[Tl are all ellipses 
whose foci are located at precisely the positions of the turning points. All of the classical orbits 
are periodic and all orbits have the same period T = 2tt. Thus, a classical particle travels faster 
along more and more distant ellipses. 
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Figure 1. Classical trajectories in the complex-x plane for the harmonic oscillator whose Hamiltonian 
\s H = + . These trajectories represent the possible paths of a particle whose energy is i? = 1. The 
trajectories are nested ellipses with foci located at the turning points a,t x = ±1. The real line segment 
(degenerate ellipse) connecting the turning points is the usual periodic classical solution to the harmonic 
oscillator. All closed paths have the same period 2-k. 

Now suppose that a classical particle of energy E = 1 \s traveling along the real-x axis from 
some point x = —a to x = a, where a > 1. If the potential V is everywhere zero along its path, 
then it will travel at a constant velocity. However, if the particle suddenly finds itself in the 
parabolic potential V{x) = just as it reaches the turning point at x = — 1 and it suddenly 
escapes the influence of this potential at x = 1, then the time of flight from x = — atox = a 
will be changed because the particle is not traveling at constant velocity between the turning 
points. Now imagine that the potential T^(x) = x^ is suddenly turned on before the particle 
reaches the turning point at x = — 1. In this case, the particle will follow one of the elliptical 
paths in the complex plane around to the positive real axis. Just as the particle reaches the 
positive real axis the potential is turned off, so the particle proceeds onward along the real axis 
until it reaches x = a. This trip will take less time because the particle travels faster along the 
ellipse in the complex plane. 

We have arrived at the surprising conclusion that if the classical particle enters the parabolic 
potential Vix) = x^ immediately after it begins its voyage up the real axis, its time of flight will 
be exactly half a period, or vr. Indeed, by traveling in the complex plane, a particle of energy 
E = 1 can go from the point x = —a to the point x = a in time vr, no matter how large a is. 
Evidently, if a particle is allowed to follow complex classical trajectories, then it is possible to 
make a drastic reduction in its time of flight between two given real points. 

2 Quantum brachistochrone problem 

The purpose of this paper is to show that by using complex non-Hermitian Hamiltonians, we 
can obtain a faster time of flight than is possible with Hermitian Hamiltonians. The quantum 
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brachistochrone problem is defined as follows: If we are given an initial quantum state |/) 
and a final quantum state \F), then there exist many Hamiltonians H under which \I) evolves 
into \F). The quantum brachistochrone problem is to find the particular Hamiltonian H that 
achieves this transformation in the least time r, subject to the constraint that the difference 
between the largest and smallest eigenvalues of H is held fixed. For Hermitian Hamiltonians, 
r has a nonzero lower bound. However, we will see that for non-Hermitian "PT-symmetric 
Hamiltonians satisfying the same energy constraint, r can be made arbitrarily small. 

One might think that this result could violate the time-energy uncertainty principle. However, 
we will see that this is not the case because for non-Hermitian Hamiltonians the path from |/) 
to l-F) can be made arbitrarily short. The mechanism is similar to that in general relativity, 
where the distance between two space-time points can be made small if they are connected by 
a wormhole. 

3 Review of quantum mechanics 

Based on the traditional training that one receives in a quantum mechanics course, one would 
expect a theory defined by a non-Hermitian Hamiltonian to be unphysical because the energy 
levels would most likely be complex and the time evolution would most likely be nonunitary (not 
probability-conserving). However, theories defined by a special class of non-Hermitian Hamil- 
tonians called T^T-symmetric Hamiltonians can have positive real energy levels and can exhibit 
unitary time evolution. Such theories are acceptable quantum theories. In principle, these theo- 
ries can be distinguished experimentally from those defined by Hermitian Hamiltonians because 
non-Hermitian time evolution can proceed arbitrarily rapidly. 

We use the following notation in this paper: By the term Hermitian, we mean Dirac Her- 
mitian, where the Dirac Hermitian adjoint symbol ^ represents combined matrix transposition 
and complex conjugation. The parity operator V performs spatial reflection x — > —x and 
the antilinear time-reversal operator T performs combined time reversal t ^ —t and complex 
conjugation. 

The first PT-symmetric quantum-mechanical Hamiltonians were introduced in 1998 [2]; and 
since then there have been many follow-up papers by a wide range of authors. There have also 
been three recent review articles O IH [5] . In [2] it was discovered that even if a Hamiltonian 
is not Hermitian, its energy levels can be all real and positive so long as the eigenfunctions are 
symmetric under VT reflection. 

These new kinds of Hamiltonians are obtained by deforming ordinary Hermitian Hamiltonians 
into the complex domain. The original class of "PT-symmetric Hamiltonians that was proposed 
in [2] has the form 



where e is a real deformation parameter. Two particularly interesting special cases are obtained 
by setting e = 1 to obtain H = + ix^ and by setting e = 2 to obtain H = — x^. Surprisingly, 
these Hamiltonians have real, positive, discrete energy levels even though the potential for e = 1 
is imaginary and the potential for e = 2 is upside-down. The first proof of spectral reality and 
positivity for H va. ^ was given by Dorey et al. in [6]. 

The philosophical background of VT quantum mechanics is simply this: The axiom of quan- 
tum mechanics that requires the Hamiltonian H to be Dirac Hermitian is distinct from all of the 
other axioms because it is mathematical in character rather than physical. The other axioms are 
stated in physical terms; these other axioms require locality, causality, stability and uniqueness 
of the vacuum state, conservation of probability, Lorentz invariance, and so on. The condition 
of Dirac Hermiticity H = is mathematical, but the condition of VT symmetry H = 
(space-time reflection symmetry) is physical because V and T are elements of the Lorentz group. 
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The spectrum in ([T]) is real, which poses the question of whether this Hamiltonian specifies 
a quantum-mechanical theory. That is, is the theory specified by H associated with a Hilbert 
space endowed with a positive inner product and does H specify unitary (norm-preserving) time 
evolution? The answer to these questions is yes. Positivity of the inner product and unitary 
time evolution was established in [7j for quantum-mechanical systems having an unbroken VT 
symmetry and in [8] for quantum field theory. 

To demonstrate that the theory specified by the in dT} is a quantum-mechanical theory, 
we construct a linear operator C that satisfies the three simultaneous algebraic equations [7]: 

= 1, [C,VT] = 0, and [C,-ff] = 0. Using C, which in quantum field theory is a Lorentz 
scalar [9], we can then construct the appropriate inner product for a "PT-symmetric Hamiltonian: 
{a\b) = aP'^'^ ■ b. This inner product, which uses the CVT adjoint, has a strictly positive 
norm: {a\a) > 0. Because H commutes with both VT and C, H is self-adjoint with respect 
to CVT conjugation. Also, the time-evolution operator e~^^^ is unitary with respect to CVT 
conjugation. Note that the Hilbert space and the CVT inner product is dynamically determined 
by the Hamiltonian itself. 

We have explained why a "PT-symmetric Hamiltonian gives rise to a unitary theory, but 
in doing so we raise the question of whether "PT-symmetric Hamiltonians are useful. The 
answer to this question is simply that PT-symmetric Hamiltonians have already been useful 
in many areas of physics. For example, in 1959 Wu showed that the ground state of a Bose 
system of hard spheres is described by a non-Hermitian Hamiltonian [10]. Wu found that the 
ground-state energy of this system is real and he conjectured that all of the energy levels were 
real. Hollowood showed that the non-Hermitian Hamiltonian for a complex Toda lattice has 
real energy levels [Tl]. Cubic non-Hermitian Hamiltonians of the form H = A- ix^ (and also 
cubic quantum field theories having an imaginary self-coupling term) arise in studies of the 
Lee~Yang edge singularity [12] and in various Reggeon field-theory models |13j . In all of these 
cases a non-Hermitian Hamiltonian having a real spectrum appeared mysterious at the time, but 
now the explanation is simple: In every case the non-Hermitian Hamiltonian is VT symmetric. 
Hamiltonians having VT symmetry have also been used to describe magnetohydrodynamic sys- 
tems |14] and to study nondissipative time-dependent systems interacting with electromagnetic 
fields [E]. 

An important application of VT quantum mechanics is in the revitalization of theories that 
have been abandoned because they appear to have ghosts. Ghosts are states having negative 
norm. We have explained above that in order to construct the quantum-mechanical theory 
defined by a PT-symmetric Hamiltonian, we must construct the appropriate adjoint from the 
C operator. Having constructed the CVT adjoint, one may find that the so-called ghost state 
is actually not a ghost at all because when its norm is calculated using the correct definition of 
the adjoint, the norm turns out to be positive. This is precisely what happens in the case of the 
Lee model. 

The Lee model was proposed in 1954 as a quantum field theory in which mass, wave- function, 
and charge renormalization could be performed exactly and in closed form [16]. However, in 
1955 Kallen and Pauli showed that when the renormalized coupling constant is larger than 
a critical value, the Hamiltonian becomes non-Hermitian (in the Dirac sense) and a ghost state 
appears [17] . The appearance of the ghost was assumed to be a fundamental defect of the Lee 
model. However, the non-Hermitian Lee-model Hamiltonian is VT symmetric and when the 
norms of the states of this model are determined using the C operator, which can be calculated 
in closed form, the ghost state is seen to be an ordinary physical state having positive norm |18j . 
Thus, the following words by Barton [19] are not true: "A non-Hermitian Hamiltonian is unac- 
ceptable partly because it may lead to complex energy eigenvalues, but chiefly because it implies 
a non-unitary S matrix, which fails to conserve probability and makes a hash of the physical 
interpretation." 
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Another example of a quantum model that was thought to have ghost states, but in fact does 
not, is the Pais-Uhlenbeck oscillator model [20]. This model has a fourth-order field equation, 
and for the past several decades it was thought (incorrectly) that all such higher-order field 
equations lead inevitably to ghosts. Indeed, it is explained in [20] when the Pais-Uhlenbeck 
model is quantized using the methods of VT quantum mechanics, it does not have any ghost 
states at all. 

There are many potential applications for VT quantum mechanics in areas such as particle 
physics, cosmology, gravitation, quantum field theory, and solid-state physics. These applica- 
tions are discussed in detail in the recent review article [1]. 

Having established the validity and potential usefulness of VT quantum mechanics, one may 
ask why VT quantum mechanics works. The reason is that CV is a positive operator, and 
thus it can be written as the exponential of another operator Q: CV = e®. The square root 
of can then be used to construct a new Hamiltonian H via a similarity transformation 
on the "PT-symmetric Hamiltonian H: H = e~2/^ife®/^. The new Hamiltonian H has the 
same energy eigenvalues as the original Hamiltonian H because a similarity transformation is 
isospectral. Moreover, H is Dirac Hermitian [21]. VT quantum mechanics works because there 
is an equivalence between a non-Hermitian PT-symmetric Hamiltonian and a conventional Dirac 
Hermitian Hamiltonian. 

There are a number of elementary examples of this equivalence, but a nontrivial illustration 
is provided by the Hamiltonian in ([T|) at e = 2, which is not Hermitian because boundary 
conditions that violate the L?' norm must be imposed in Stokes wedges in the complex plane in 
order to obtain a real, positive, discrete spectrum. The exact equivalent Hermitian Hamiltonian 
\s H = + — 2hx^ where h is Planck's constant [22]. The term proportional to h vanishes 
in the classical limit and is thus an example of a quantum anomaly. 

We have established that VT symmetry is equivalent by means of a similarity transformation 
to conventional Dirac Hermiticity. Therefore, one may wonder whether VT quantum mechanics 
is actually fundamentally different from ordinary quantum mechanics? The answer is yes, and 
this paper argues that, at least in principle, there is an experimentally observable difference 
between PT-symmetric and ordinary Dirac Hermitian Hamiltonians. The quantum brachis- 
tochrone provides a setting for examining this difference and provides a way to discriminate 
between the class of PT-symmetric Hamiltonians and the class of Dirac Hermitian Hamilto- 
nians. 

4 Solving the Hermitian quantum brachistochrone problem 

To find the Hermitian Hamiltonian H that solves the quantum brachistochrone problem we 
must examine all possible Hamiltonians under which a state in Hilbert space evolves into 
another state \'iPf) in time t: 

The problem is to find the minimum time t = t required for this transformation, subject to the 
constraint that the difference lo between the largest and smallest eigenvalues of i7, 

is held fixed. The quantum brachistochrone is the Hamiltonian that performs this time evolution 
in the least possible time. In [23] it is shown that for Hermitian Hamiltonians r 7^ 0. However, we 
show in this paper that one can find a Hamiltonian in the space of PT-symmetric Hamiltonians 
that satisfies the same energy constraint and can perform the time evolution in no time at all! 
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Here, we study the simplest case of Hamiltonians having only two energy levels. We restrict 
the discussion to this case because it is shown in [21] that one need only work in the two- 
dimensional subspace of the full Hilbert space that is spanned by the initial state vector l-i/^/) 
and the final state vector \iPf)- We consider the case of Hermitian Hamiltonians and choose 
a basis so that 

m = ( J ) and = ( I 

where the condition that \iPf) be normalized is |ap + = 1. The most general 2x2 Hermitian 
Hamiltonian is 

H = { ^jfj ) (r, s, u, 9 real), 

y re u J 

For this Hamiltonian the eigenvalue constraint takes the form 

= (s - n)2 + 4r2. (2) 

To find the optimal Hamiltonian satisfying this constraint, we express H in terms of the 
Pauli matrices: 

H = |(s + u)l + iwfT-n, 

where 

n = — (2r cos 9, 2r sin 9, s — u) 

is a unit vector and 

01\ f Q -i\ /lO 



^1 = ' 1 ' ^2 = U ' V -1 



We use the identity 

exp(i(/)(T-n) = cos (/) 1 + i sin (/) (T • n 
to write IV'f) = e~^^'^/'^\ipi) as 



\ UJ Zfl / 

The second component of this equation then gives |6| = ^sin^, which allows us to find the 
required time of evolution: 

2h uj\b\ 

t = — arcsin . 

oj 2r 



We must now minimize the time t over all r > while maintaining the constraint in ([2]) . This 
nstraint tells us that the maximum value 
minimum evolution time r is thus given by 



constraint tells us that the maximum value of r is ^u). At this maximum we have s = u. The 



TOJ = 2/i arcsin \ h\. 



(3) 
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Note that if a = and b = 1, we have r = nh/uj for the smallest time required to transform 
( Q ) to the orthogonal state ( ^ ) . The time r required to transform a vector into an orthogonal 
vector is called the passage time. 

The form of the result in ([3|) resembles the uncertainty principle, but ([3]) is merely the 
statement that rate x time = distance. The constraint in ([2]) is equivalent to a bound on the 
standard deviation AH, where (AH)'^ = — {Tp\H\'ip)'^ in the normalized state lip). The 

maximum of AH is lv/2. The speed of evolution of a quantum state is given by AH. The 
distance between the initial state \tp[) and the final state IV'f) is 2 aiccos{\{ipF\tpi)\) ■ Thus, the 
time T to evolve from \^pJ) to {iPf) is bounded below because the speed is bounded above with 
the distance held fixed. 



5 Solving the non- Hermitian quantum brachistochrone problem 

For a "PT-symmetric Hamiltonian, r can be arbitrarily small. This is because a PT-symmetric 
Hamiltonian whose eigenvalues are all real is equivalent to a Hermitian Hamiltonian via H = 
e-2/2i/e2/2. The states in a T^T-symmetric theory are mapped by e to the corresponding 
states in the Dirac Hermitian theory. But, the overlap distance between two states does not 
remain constant under a similarity transformation. We can exploit this property of the similarity 
transformation to overcome the Hermitian lower limit on the time r. The detailed calculation 
is explained in [25], and this calculation has already led to much research activity and lively 
debate [26]. 

We consider the general class of PT-symmetric 2x2 Hamiltonians having the form 

^ = ( 'f J-^e ) ir, s, e real), (4) 
where T is complex conjugation and V = q). The eigenvalues 



E± = r cos 6 lb \/ s"^ — r"^ sin^ 9 

are real if > sin^ 6. This inequality defines the region of unbroken VT symmetry. The 
unnormalized eigenstates of H are 

oia/2 \ / A„-ia/2 



where a (real) is given by sin a = (r/s) sin^. The C operator for H m ^ is 
1 / zsina 1 



cos a 



1 —2 sin a 



It is easy to verify that the CVT norms of both eigenstates have the value \/2 cos a. 
To calculate r we express the H in ([4]) as 

H = {r cos 9)1 + ^uja ■ n, 

where 

2 

n = — (s, 0, ir sin 9) 
is a unit vector. The squared difference between energy eigenvalues is 



2 



As"^ -4r'^sin'^9. (5) 



The positivity of a;^ is ensured by the condition of unbroken VT symmetry. 
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To determine r we write down the T^T-symmetric time-evolution equation: 

V y cosa -«sin(g) ) ' 

Consider the pair of vectors used in the Hermitian case: \'>^i) = (q) and \iIjf) = (?)• (Note that 
these two vectors are not orthogonal with respect to the CVT inner product.) Observe that the 
evolution time needed to reach \iPf) from l^/^/) is t = {2a — ■K)h/uj. Optimizing this result over 
allowable values for a as a approaches ^vr, the optimal time r tends to zero! 

6 Discussion 

Equations ([2]) and ([5]) reveal the difference between Hermitian and PT-symmetric Hamiltoni- 
ans. Equation ([2]) for the Hermitian matrix Hamiltonian has a sum of squares while ([5]) has 
a difference of squares. The elliptic equation ([2]) gives a nonzero lower bound for r. The hyper- 
bolic equation dS]) allows r to approach zero because the matrix elements of a "PT-symmetric 
Hamiltonian can be made large without violating the energy constraint Ej^ — = uj. The fact 
that T can be made arbitrarily small may have applications in quantum computing. 

We conclude with two comments. First, as a — > ^vr we get cosa — > 0. However, the energy 
constraint becomes up' = 4s^ cos^ a. Since uj is fixed, to have a approach ^vr, we must require 
s ^ 1. It follows from sin a = (r/s)sm6 that |r| ~ so we must also require that r ^ 1. 
Thus, if T <C 1, the matrix elements of the "PT-symmetric Hamiltonian are large. Second, the 
result that r = does not violate the uncertainty principle. Both Hermitian and non-Hermitian 
'PT-symmetric Hamiltonians share the properties that (i) the passage time is given by nh/u;, 
and (ii) AH < io/2. 

To summarize, the key difference between the Hermitian and the non-Hermitian case is that 
(o) and (?) are orthogonal in the Hermitian case, but they have separation vr — 2|a| in the 
PT-symmetric case. This is because the Hilbert space metric of a PT quantum theory depends 
on H. By choosing the parameter a properly, we create a wormhole-like effect in Hilbert space. 
That is, we find a path in Hilbert space from the initial state vector It/jj) to the final state 
vector \iPf) that is shorter than the Hermitian path. This is analogous to finding a wormhole 
in coordinate space. In short, what we have done here is to construct a "wormhole" in Hilbert 
space. 
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